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1. INTRODUCTION 


Clouds play a significant role in remote sensing of the earth's features from 
altitude by electromagnetic means. Depending on its composition and thickness, a 
cloud can partially attenuate or totally block the passage of the electromagnetic radi- 
ation from the target to an airborne or spaceborne sensor. With judicious choice of 
the wavelength of the radiation to be measured, these attenuating effects may be 
minimized, but never eliminated. However, the nature of the experimental measure- 
ment itself may dictate the selection of a particular spectral interval such that no 
meaningful tradeoffs can be made between target identification and cloud attenuation. 
Under these circumstances, knowledge of the global spatial and temporal distributions 
of cloud parameters is essential in the successful design of remote sensors and 
remote sensing experiments. In recognition of this fact, the Aero-Astrodynamics 
Laboratory initiated a program to establish a statistical data bank of cloud parameter! 
on a global basis, and to develop computer techniques whereby these statistics may 
be used for cloud cover simulation experiments. The results presented in this 
report represent the laieui liiiprOvemerits in the cloud statist^ data bank and in the 
simulation scheme, 

1. 1 Objectives 

The basic objective of the study reported in this document was the improve- 
ment of the existing cloud data bank by the addition of statistics of cloud type, 
cloud layers and the microphysical properties of clouds. More importantly, these 
statistics were to be formatted such that they could be incorporated into an improved 
simulation scheme to provide not just a statement of the probability of encountering 
a given cloud cover condition at a specified time and place, but also the probable 
cloud type, drop size, liquid water content and other parameters which ultimately 
determine the degree of electromagnetic attenuation. 

A second objective of the study was to develop a means of adjusting the 
current simulation output as a function of the ratio of cloud size to sensor reso- 
lution size. This objective, in detail, appears unrelated to the more important 
tasks of developing additional statistical parameters to characterize the cloud 
field. However, its achievement does represent an overall improvement in the 
whole simulation program. 
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1.2 specific Tasks 


In meeting the overall objectives discussed previously* four specific tasks 
to be performed were defined at the outset of the study. For completeness, these 
tasks are itemized below: 

Task 1 

Develop a means of adjusting the current simulation output as a 
function of the ratio of cloud size to sensor resolution size. 

Task 2 

Develop cloud -type frequency distributions by month for each 
homogeneous cloud region. 

Task 3 

Develop a statistical model to correlate cloud type with cloud 
thickness, water content, and water drop size distribution. 

Task 4 

Develop cloud thickness, water content, and water drop size 
statistics for each homogeneous region. 

It is clear from the above that the last three tasks, taken together, are intended 
to establish the global statistical distributions of microphysical cloud parameters. 

1.3 Summary of the Contents of the Report 

In the study, primary effort was directed towards Task 2, This was made 
necessary by the fact that while a great deal of raw data exist on cloud type obser- 
vations, very little has been done in the compilation of detailed climatological sum- 
maries. Since cloud-type and cloud-layer statistics form the foundation upon which 
the whole data bank of global distributions of microphysical cloud properties is based, 
it became an important and necessary task to process existing cloud observations, 
in the form of WBAN entries, to provide these essential statistics. The procedures 
to derive these statistics, and the statistics themselves, occupy, therefore, a large 
portion of this report. 
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The statistical results reported in this document have significance in them- 
selves, However, the format in which they are presented has greater significance 
in the context of the specific problem of cloud simulation for remote sensing obser- 
vations, For this reason, a short recapitulation is provided in the next section of 
past studies in cloud simulation initiated and directed by the Aero -Astr ©dynamics 
Laboratory. The following Sections 3 and 4 discuss the sources of the data and the 
procedures used in the current study to compile the cloud statistics. Some of the 
problems in deriving the statistics are also discussed. Examples of the output data 
formats for recording the statistics are shown and documented. 

In Section 5, the results of cloud-type statistics are discussed. This discus- 
sion is directed towards those features in the statistics of particular importance to 
the use of the data in cloud simulation programs, including such topics as regional 
homogeneity, sample size dependence, temporal variability and regional dependence. 
The results and format of the statistical compilation of the microphysical properties 
of clouds are separately discussed in Section 6, In Section 7, the work on develop- 
ment of cloud statistics is concluded with a discussion of how the new statistics may 
be com bine a with the existing data bank in a multi step Monte Carlo simulation schema 


for earth surveillance. 

The results of Task 1 on resolution are separately presented as an appendix. 
This organization was arrived at so as not to break the continuity of the discussion 
of the cloud -type statistics. The actual statistics of cloud type, layer and micro - 
physical properties compiled by the study are tabulated in a second appendix. 
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Z, BACKGROUND 


The basis of the current study is contained in two previous studies initiated 
by the Aero-Astrodynamics Laboratory. Since the results presented in this report 
can best be understood in the context of the past studies, summary discussions of 
their results are presented below. 

2, I Worldwide Cloud Cover Distributions for Use in Computer Simulation 

This first study by Sherr et al (1968) defined the basic purpose for which 
cloud statistics are compiled and laid down a number of ground rules for the types 
and format of the statistics themselves. These ground rules were further reinforced 
by the particular simulation scheme developed for the application of the statistics. 

As defined by Sherr et al (1968), the primary purpose for the development 
of a cloud statistics data bank was to provide data for proper mission analysis and 
simulation during the planning of future missions, and to determine the probable 
success of already planned earth-oriented space missions. In this light, they 
considered that simulation requires a set of cloud climatological data of fairly 
universal application. While uncertainties in simulation preclude the necessity of 
extreme accuracy, the data bank was to provide at least the following: 

1) Global coverage 

2) Cloud cover distributions in a readily usable, standard form. 

3) Distributions by season, time of day, and some readily defined 
climatological region or grid. 

4) Expression of the spatial and temporal coherence of cloud cover. 

5) The expression of cloud cover distributions on a variety of scales 
of observation. 

To achieve these design objectives, statistical analysis and compilation were per- 
formed of cloud amounts for monthly, seasonal and annual reference periods, for 
various times of day, for selected stations which typify the various cloud climatolo- 
gies over the globe. In addition to these unconditional cloud cover probabilities, 
temporal and spatial conditional cloud cover statistics were also compiled for the 
selected stations. 
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Table 2-1 


Representative Stations for Unconditional Distributionc 


REGION 

NUMBER 

STATION 

COORDINATES 

TYPE OF 
DATA 

YEARS oy 
RECORD 


J 

Dhahran. Saudi Arabia 

26-17N 50-09E 

2 

11 



(Airfield) 






Tripoli, Libya 

32-54N 13-17E 

1 

19 



(Wheelus AFB) 





i 

Angeles, lAiSon, P. I. 

15-llN 120-33E 

1 

21 



(Clark AFB) 





•4 

Tampa, Florida 

27-51N 82-30W 

1 

23 



(M.icDill AFB) 





1 

5 

L Angeles, California 

33-56N 118-23W j 

1 

19 


1 

(WHAS) Hours 10-19 (May-October) 

i 


j 

i 

Modified 



1 

i 

Talara, Peru 

04-32S B1-14W 

2 

S j 

\ 

01 and 22 Hours -Synthetic 

i 







i ’ 




i 

1 8 

Mountaix. Home, Idaho 

43-03N 115-51W 

1 

20 ! 


(AFB) 




! 9 

Fort Yukon, Alaska 

66-35N 145-ieW 

2 

18 i 


(WB) 




' 10 

Harbin, China 

45-45N 126-38W 

3 

7 

11 

Belleville, Illinois 

38-33N 89-51W 

1 

- ' 


(Scott AFB) 



. 

12 

Ban Me Thuot, Vietnam 

12-41N 108-07E 

1 


• 

(City Airport) 




13 

Ship D 

44-OON 41-OOW 

4 

10 


(Atlantic) 



1 

M 

Adak, Alaska 

51-S3N 176-38W 

1 

2S 


(NS) 



1 ; 

IS 

Resolute, NWT, Canada 

74-41N 94-55W 

2 

1 ‘ 

i 

Fort Kobbe, Canal Zone 

08-55N 79-36W 

1 

! 19 

1 

(Howard AFB) 



f 

17 

Bangalore. India 

12-57N 77-38E 

3 

7 1 


(Hindustan Airport) 



i 

1 

San Francisco, California 

37-37N 122-23W 

1 

18 j 


(WBAS) 




19 

Shreveport, Louisiana 

3Z-30N 93-40W 

1 

27 


(Barksdale, AFB) 



1 

20 

Ship V 

31-OON 164-OOE 

4 

10 1 


(Pacific) 







1 Lesend for Type of DaU: 



21 

Seasonal Reversal of Region 12 

(1) Revised Uniform Summary (A-F) 


22 

Seasonal Reversal of Region 13 

(2) Original Uniform Summary (A and B) 


IJ 

Seasonal Reversal of Region 14 

(3) NIS Summary 





(4) Raw OaU (Ship Sutions) 



24 

Seasonal Reversal of Region 15 





25 

Seasonal Reversal of Region 16; Hours 10, 13. 16 for May - 

September Modified 

26 

Seasonal Reversal of Region 17 



j 

27 

Seasonal Reversal of Region 18 




28 

Seasonal Reversal of Region 19 




29 

Seasonal Reversal of Region 20 
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Basic to the development of the cloud statistics are the cloud climatological 
regions into which the global surface was divided. These regions are shown in 
Figure 2-1. They were designed primarily on the basis of climatological consider- 
ations. A detailed description of each of these regions may be found in the report 
by Sherr et al (1968). In addition, a representative, or prototype, station was 
selected for each of the regions. The statistics derived using data from a repre- 
sentative station were assumed to apply to the region as a whole. Table 2-1 shows 
a tabulation of the representative stations, their geographical locations, and the 
type and length of record used to generate the unconditional cloud cover statistics. 

Cloud cover situations were considered to fall into only the five categories 
shown in Table 2-2. The corresponding cloud cover specifications in tenths and in 
octas as commonly used in reporting cloud cover are also shown. Using the cloud 
summaries indicated in Table 2-1, cloud cover statistics for each of the five cate- 
gories were derived for every one of the regions, for twelve months of the year, and 
at three-hour intervals throughout the 24-hour day beginning with 01 LST. Spatial 
and temporal conditional statistics of cloud cover were also derived using primarily 
meteorological satellite photographs. The results of these statistical compilations 
were preserved on punched cards for computer simulation procedures. 


TABLE 2-2 

CLOUD CATEGORY DESIGNATION 


Category 

Tenths 

Eighths 
(Octas ) 

1 

0 

0 

2 

1,2,3 

1, 2 

3 

4. 5 

3,4 

4 

6. 7, 8, 9 

5,6,7 

5 

10 

8 


c 
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As a final task in the studyt a Monte Carlo technique was developed whereby 
the unconditional and conditional statistics, in the form of cumulative probabilities, 
are used to simulate the cloud cover conditions for various orbital missions. With 
some modifications this simulation technique has remadned the basic approach 
towards the application of the statistical data. 

2.2 Development of a Global Cloud Model for Simulating Earth- Viewing 
Space Missions 

This second study, by Greaves et al (1971) continued the development of 
statistics and statistical procedures for cloud simulation in a number of areas. 

Most important of these is the development of a Markov scaling technique to scale 
the 24-hour and 200-mile conditional cloud cover statistics to any other time or 
distance scale. Also developed was a relationship between ground and satellite 
derived cloud frequency distributions. This relationship was used to introduce 
internal consistency between the ground-observed unconditional statistics and the 
conditional statistics derived from -satellite observations. 

Also in this study, a first attempt was made to establish some statistics 
for the occurrence of weather phenomena, other than clouds, such as tornados and 
severe thunderstorms which are of significance to remote sensing of the earth* s 
surfaces from altitude. 

2.3 The Role of the Current Study 

It is clear, even from the limited summaries provided above, the previous 
studies concentrated only on the cloud cover aspects of the total cloud attenuation 
problem. For sensors such as cameras and photometers operating in the visible 
part of the electromagnetic spectrum, the specification of percentage of cloud 
cover is adequate to describe the expected observational coverage. However, for 
sensors operating in the longer wavelengths, such as millimeters and centimeters, 
clouds need not be totally interfering. A more detailed specification of the cloud 
field is required for the determination of attenuating effects. Parameters such as 
cloud thickness, drop size distribution, liquid water content and number density 
become important. It is towards the development of these statistical parameters 
that the current study is addressed. 
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3, DATA SOURCE 


3.1 Data Search 

A systematic search was conducted to locate existing cloud summaries or 
cloud climatologies which might provide useful cloud-type and cloud-layer statistics 
for one or more of the cloud regions. In this search the aid of the personnel at the 
National Climatic Center (NCC) was sought and received. They identified a total 
of nineteen studies on file at the Center with data related to the distribution of cloud 
types. The su m mary file cards for these studies v/ere obtained and examined. An 
example of the information contained in the file cards is shown in Table 3-1. It is 
clear from this typical file card that the summary data on record with the NCC are 
too specialized and not appropriate for the task at hand. 

The search conducted by ARA personnel located a number of sources which 
included the following examples: 

Crutcher, H. L. and O. M. Davis. 1969: U. S. Navy Marine Climatic 

Atlas of the World, Vol. VII: The World, NAVAIR 50-1C-54, by 
Direction of Commander, Naval Weather Service Command. 

In this document, marine cloudiness is represented by a series of three 
charts: Total Cloudiness, Low Clouds, and Median Cloudiness, derived from 

analysis of 1200 GMT observations made at ship and coastal stations. 

McDonald, W, F. , 1938: Atlas of Climatic Charts of Oceans, W. B. 

No. 1247, Weather Bureau. U. S. Dept, of Agriculture. 

In this atlas, average frequencies of occurrence by season of six cloud types 
from ocean regions between SO^'S to 6 o"n are provided. The cloud-type categories 
are: (1) Ci. Cs, Cc; (2) As. Ac; (3) Cu; (4) St. Sc; (5) Cb; and (6) Ns. 

On a more local scale we located a series of reports prepared by the Institute 
of Atmospheric Physics, the University of Arizona, which analyzed in detail cloud 
field parameters over Arizona. Typical of these are: 

Des Jardins. R. B. . 1958; The Distribution of Cl ouds at Tucson, 

Arizona, with Respect to Type, Amount and Time of Observation, 
Scientific Report No. 6 under Contract NSF G- 1101, Institute of 
Atmospheric Physics, University of Arizona. 

Sellers, W. D, , 1958: The Annual and Diurnal Variations of Cloud 

Amounts and Cloud Types at Six Arizona Cities, Scientific Report 
No. 8, University of Arizona, Institute of Atmospheric Physics. 
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TABLE 3-1 


EXAMPLE OF CLOUD SUMMARIES ON FILE AT NCC 


Card Serial Number; A-372 

Remarks: This is percentage frequency tabulation of specified 

cloud conditions by st-mo as shown: 


where 


1. Sky clear, or A except when foredorninant type I 

2. Total elds C when predominant type V 

3. Total elds C when predominant type IV 

4. Total elds C when predominant type III 

5. Total elds C when predominant type II, spread by 5000 ft 

heights . 

6. Total elds A when predominant type is I 

7. Total elds B when predominant type is I 

8. Total elds B when predominant type is V 

9. Total elds B when predominant type is IV 

10, Total elds P when predominant type is III 

11, Total elds B when predonninant type is II 

Elements 
Cloud - amount 

Cloud - type 


Class Intervals 

A, 0-3 tenths 

B, 3-6 tenths 

C, 6-10 tenths 

I Cb, Cu, Sc, Ace, Cm 

II St, Fs, Fc, Ns 

III Cs, Ac 

IV Cs, Cc 

V Ci 

VI Obscurations 
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Of the existing cloud summaries, none was found totally satisfactory. The 
detailed summaries of Arizona compiled by the University of Arizona perhaps come 
closest to providing the data to generate the necessary statistics, at least for Region 
2. Even these data are deficient in cloud layer statistics. The Atlas of Climatic 
Charts of the Oceans provides the most complete areal coverage of the ocean areas. 
However, its deficiencies are self-e\ddent in that the data are compiled solely from 
1200 GMT observations and cover only the ocean areas with no detailed breakdown of 
the cloud types. Thus for the vast Pacific Ocean areas, where a large diurnal effect 
is known to take place, cloud-type statistics are available only during the nighttime 
hours where cumuliform clouds are at an unrepresentative minimum, at least unrep- 
resentative of the true daily mean statistics. 

In the light of the inadequacies of existing summaries, it was felt that the 
only way to obtain these statistics was to generate them from the raw surface 
observations. In a study of the cloud climatology off the southern California coast, 
Greaves et al (1969) had in fact developed limited cloud-type statistics by process- 
ing over 10 years of WBAN data recorded on tape. The experience gained in this 
study showed that the WBAN records on tape provide the ideal raw data for the 
development of cloud -type and cloud -layer statistics. 

3,2 Selection of Representative Stations 

The global statistics of cloud types and layers are composited from statistics 
derived for each of the 29 cloud climatological regions into which the surface of the 
globe was divided. These regions are exactly those previously developed and 
defined by Sherr et al (1968) in their derivation of global statistics of cloud cover. 
Since in any cloud simulation scheme, cloud cover itself still remains the para- 
mount cloud parameter, redefinition of the cloud climatological regions on the 
basis of cloud-type or cloud-layer statistics was not considered necessary. 
Furthermore, it was the intent to make use of the same set of representative 
stations, used by Sherr et al, to develop the new statistics. By so doing, self- 
consistent sets of cloud statistics would be obtained. A search was therefore 
initiated to locate and obtain taped WBAN records for each of the representative 
stations previously identified in Table 2-1. 

In the search for long period WBAN records, it was ascertained from the 
National Climatic Center (NCC) that taped observations were available for the 
following cloud regions. Some of these tapes were already on file at ARA. 
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TABLE 3-2 


TAPED DATA STATIONS 


Cloud Region 

Station 

1 

Dhahran, Saudi-Arabia 

2 

Tripoli, Libya (Wheelus AFB) 

3 

Angeles, Luzon, P. I. (Clark AFB) 

4 

Tampa, Florida (MacDill AFB) 

5 

Los Angeles, California (WBAS) 

8 

Mountain Home, Idaho (AFB) 

9 

Fort Yukon, Alaska (FAA) 

11 

Belleville, Illinois (Scott AFB) 

13 

Ocean Station D 

14 

Adak, Alaska 

16 

Fort Kobble, Canal Zone (Howard AFB) 

18 

San Francisco, California (WBAS) 

19 

Shreveport, Louisiana (Bardsdale AFB) 


20 


Ocean Station V 




It is clear from Table 3-2 that taped WBAN data are available from NCC only for 
stations operated by the United States, either by the USAF, the USWB or, in the 
case of Fort Yukon, by the FAA. 

This latter station, Fort Yukon, presented a unique problem since it was 
subsequently found that the tape from Fort Yukon in fact contained no cloud -type 
data. It was therefore necessary to select Fairbanks, Alaska, as the new repre- 
sentative station for Region 9. Taped data were available for this station. No 
taped data were located for representative stations in each of the remaining six 
stations (Regions 21 to 29 being seasonal reversals of Regions 1 to 20). For these 
regions, WBAN data are available, if at all, only as microfilm records. In con- 
sultation with the staff at the NCC, appropriate stations were selected to provide 
representative statistics for the remaining regions. These are shown in Table 3-3. 

TABLE 3-3 

MICROFILM RECORD STATIONS 


— 

Region 

Station 

6 

Talara, Peru 

10 

Felly Bay, N, W. T. 

12 

Ton Son Nhut, Vietnam 

15 ' 

Resolute, N. W. T. 

17 

Bangalore, India 


It was impossible to locate station data for Region 7 which was formulated originally 
by Sherr et al (1969) on the basis of synthesized data. As a first attempt at providing 
cloud-type and layer statistics for this region, it was decided to extend the contiguous 
Region 5 statistics into Region 7. This extrapolation may not be totally valid. 
However, it was noted by Sherr et al that Region 7 and Region 5 differ only in some 
seasonal aspects of the cloud cover. Until new sources of data can be identified 
for Region 7, the extrapolation adopted will have to siiffice. 

Felly Bay was selected for Region 10 instead of the original Region 10 
station Harbin, China. This was made necessary because of the unavailability of 
WBAN data from China, Ton Son Nhut replaced Ban Me Thuot as a matter of con- 
venience since data from the latter station were already on file. Three other 
changes were made in the representative stations. Because we already had on file 
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taped WBAN data from San Diego (Region 5), Ship K (Region 13) and Ship C (Region 
14), these stations were used instead of Los Angeles (Region 5), Ship D (Region 13), 
and Adak (Region 14), By so doing the time and expense needed to secure additional 
data were saved. (For a complete list of the representative stations, see Table 3-4.) 

3.3 Periods of Record 

The initial intent was to establish the cloud-type and cloud-layer statistics 
using 10 years of data. For this reason, taped WBAN records for the 10-year 
period, January 1958 to December 1967, were obtained for Regions 2, 3, 4, 5, 8, 

11, 16, 18, and 19. For Regions 1 and 9 data from the identical 10-year period 
were not available. Instead, for Region 1, the 10-year record period used spanned 
January 1952 to December 1961; and for Region 9, the period was from January 
1953 to December 1962. Subsequent to the receipt of the data, it was found that at 
USAF-operated stations, the practice of recording cloud-type data on punched cards 
was terminated as of 1 January 1964, This in effect reduces the effective length of 
the records for these stations to six years. 

At the outset of the study, it was recognized that to reduce 10 years of 
WBAN observations recorded on microfilm by manual means would be time-con- 
suming and costly even for one station. As noted in Section 3, 1. I, microfilm data 
had to be reduced for five stations. Within the time and financial limitations of the 
study, the task would be clearly impossible. It was therefore decided to provide 
only a first attempt at establishing the cloud statistics for these regions by confining 
the analysis to no more than three years' record for each of the regions. In fact, 
when the data were actually requested from the NCC, it was found that for two of 
the stations (Bangalore and Talara) the actual record length was limited to a few of 
the years of World War II between 1941 and 1945. 

Table 3-4 provides a summary of the representative stations and the lengths 
of records used in this study to extract cloud-type and cloud-layer statistics. Also 
included for comparison are the tabulations of stations and lengths of data records 
used by Sherr et al in their study of unconditional cloud cover statistics. 

3 . 4 Data F ormat 

3.4,1 Microfilm 

In the microfilm format, the cloud-type and cloud-layer data are generally 
recorded on standard WBAN-lOB forms. In this format, hourly surface observa- 
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TABLE 3-4 


REPRESENTATIVE STATIONS FOR CLOUD COVER, 
CLOUD TYPE AND LAYER STATISTICS 


Region 

Station 

Length of 
Record 
(years) 

Station 

Length of 
Record 
(years) 

1 

Dhahran (Saudi Arabia) 

11 

Dhahran (Saudi Arabia) 

6 

2 

Tripoli (Libya) 

19 

Tripoli (Libya) 

6 

3 

Angeles (Philipines) 

21 

Angeles (Philipines) 

6 

4 

Tampa, Florida 

23 

Tampa, Florida 

6 

5 

Los Angeles, California 

19 

San Diego, California 

10 

6 

Talara (Peru) 

5 

Talara (Peru) 

2 

7 

Synthetic Data 


^ 

«j»y JUbiici'j. A^ci«,o 


8 

Mountain Home, Idaho 

20 

Mountain Home, Idaho 

6 

9 

Fort Yuka, Alaska 

18 

Fairbanks, Alaska 

10 

10 

Harbin (China) 

7 

Pelly Bay (Canada) 

2 

11 

Belleville, Illinois 

j 

27 

Belleville, Illinois 

6 

12 

Ban Me Thuot (Vietnam) 

10 

Ton Son Nhut (Vietnam) 

3 

13 

Ship D 

10 

Ship K 

10 

14 

Adak, Alaska 

25 

Ship C 

10 

15 

Resolute (Canada) 

7 

Resolute (Canada) 

2 

16 

Fort Kobbe, Canal Zone 

19 

Fort Kobbe, Canal Zone 

6 

17 

Bangalore (India) 

7 

Bangalore (India) 

2 

18 

San Francisco, California 

27 

San Francisco, California 


19 

Shreveport, Louisiana 

27 

Shreveport, Louisiana 

6 

20 

Ship V 

10 

Ship V 

10 


15 








tions are recorded of sky conditions, visibility, pressure, temperature (dry and 
wet bulb), wind speed and direction, altimeter setting, and most important of all, 
detailed accounting of the type, coverage and base height of each layer of cloud 
observed. In normal practices, these observations are made every hour on the 
hour and entered into the WBAN-lOB forms by hand. However, with the data used 
in this study, this practice was not always adhered to. Furthermore, depending 
on the station, the cloud types were sometimes specified in code and not always 
recorded in WBAN forms. Some of the more important characteristics in the 
data format are tabulated in Table 3-5. This table shows clearly some of the 
problems which had to be faced in deriving cloud -type stritistics at 3 -hour intervals 
through the 24-hour day. The extrapolations used to achieve these statistics for 
stations such as Pelly Bay are discussed in Section 4, 

3.4.2 Taped Data 

Observations of cloud amount and cloud type and the manner in which they 
were formatted on magnetic tape differed somewhat between ocean vessels and 
land stations. For the land stations, a single code using 16 alphanumerics and a 
symbol is employed to describe 13 cloud types (including such types as stratus 

fractus/fractostratus and altocumulus castellanus), three obscuration phenomena, 
and the clear-sky condition. In this code, no distinction is made between high, 
middle, or low clouds. The ship reports, on the other hand, use three codes 
to separately record high, middle, and low clouds. Each code consists of 1 1 
alphanumerics (0 to 9 and X). In this scheme, the cloud type can be recorded in 
greater detail. 

For most of the stations, cloud-type data were recorded on tape only for 
the three hourly intervals corresponding to 00, 03, 06, .... GMT. As a conse- 
quence, in a number of cases it was not possible to derive the cloud-type statistics 
precisely for 01, 04, 07, .... LST. Table 3-6 s ows the actual local times of 
the cloud-type data available for the stations in q> estion. 
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TABLE 3-6 


LOCAL STANDARD TIME OF CLOUD-TYPE OBSERVATIONS 


Region 

Station 

Local Standard Time 

3 

Angeles, P. I. 

02, 05 23 

8 

Mountain Home, Idaho 

02. 05, 23 

11 

Belleville, Illinois 

DO, 03 21 

13 

Ship K 

02, 05 23 

19 

Shreveport, Louisiana 

00, 03, 21 

20 

Ship V 

02, 05 23 









4. DERIVATION OF CLOUD-TYPE AND CLOUD-LAYER STATISTICS 


4. 1 Selection of Statistical Parameters 


Given the nature of the WBAN records as described in the previous section, 

It is clear that many statistical parameters can be derived to describe the charac- 
teristics of the four -dimensional cloud field. For example, a number of conditional 
statistics can be derived. It is not unreasonable to assume that knowledge of the 
probability of occurrence of cloud type A, given that cloud type B is known to have 
occurred, is essential in certain applications. With nine basic categories of cloud 
types (Cu. Sc, St, Cb, As, Ac, Ci. Cs, Ns), at least three conditions of cloud 
layers, and the five categories of unconditional cloud cover situations, the number 
of possible combinations of the parameters into conditional statistics of the type 
suggested is extremely large. To reduce the problem to more manageable terms 
and to develop statistics which would be of direct use in cloud simulation, a simu- 
lation procedure had to be anticipated, at least in broad terms, even though this 
was not a specific task of the study. 

The procedurei as envisioned, considers the five categories of unconditional 
cloud cover statistics as basic to the simulation of the four-dimensional cloud field. 
Having specified the unconditional cloud cover category through the Monte Carlo 
technique described by Sherr et al (1968), the simulation process proceeds to specify 
the number of cloud layers, again through a Monte Carlo technique using derived 
statistics of the frequency of occurrence for 1, 2, or 3 layers of clouds, W^ith the 
percentage of cloud cover and the number of cloud layers established, the cloud 
type associated with each layer is then specified through a third application of the 
Monte Carlo technique now applied to the cloud-type statistics. The specifications 
of the three macro parameters of the cloud field are now complete. Statistical 
models of the microphysical parameters of each of the nine cloud types are then 
applied to arrive at the vertical distributions of such parameters as liquid water 
content, drop size distribution, and others which are of significance to earth obser- 
vations from altitude. 

Given the above concept of a simulation scheme, a scheme described in 
greater detail in Section 7, it is clear that the following sets of statistical values 
are needed. 
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This last set of cumulative probability values will permit the use of the Monte 
Carlo technique to determine the occurrence of multilayer clouds in simulation. 

A similar set of statistical parameters can be generated for cloud type 
as follows: 
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B.2 


Normalized Cloud-Type Frequencies or Probabilit 




where 

i 

C. = S F., and C = 1. 

‘ j=i j ^ 

Again, the last set of statistical values will permit the use of Monte Carlo techniques 
to specify cloud types in simulation procedures. 

It is important that a set each of the cloud -layer and cloud -type statistics 
be generated for each of the four nonclear-sky cloud-cover categories 2, 3, 4, and 
5, making these statistics, in fact, conditional statistics, conditioned on the cloud 
cover situation. By so doing, any dependence of the occurrence of a particular 
cloud type, or number of cloud layers, on the amount of cloud is automatically 
accounted for. Similarly, the seasonal and diurnal dependence of the cloud-type 
and cloud-layer statistics can be accounted for by obtaining the monthly statistics 
for each of the cloud regions at three -hour intervals. Thus for each of the regions, 
a total of 396 sets each of the cloud-layer and cloud -type statistics will be required. 

4.2 Data Extraction 

4.2. 1 Magnetic-Taped Data 


Two separate computer programs were developed to process the taped WBAN 
data to handle the differences in the formatting of the data between the ocean and land 
stations. These two programs were designed to accumulate frequencies of cloud types. 


21 



cloud layers and cloud cover statistics. The latter statistics, cloud cover, were 
also compiled from the tapes, even though they have already been derived by 
Sherr et al, because their compilation represented no significant increase in 
time and cost of the data -processing . Secondly, the new statistics derived from 
a different size data sample, will permit an investigation of the sensitivity of 
these statistics to sample size. More important, however, is the fact that by 
obtaining all of the cloud statistics from the same data sample, a self-consistent 
set of values will be achieved. 

In the processing, nine categories were chosen to represent the major 
cloud types. These are: cumulus (Cu), stratocumulus (Sc), stratus (St), cumulo- 

nimbus (Cb), altostratus (As), altocumulus (Ac), cirrus (Ci), cirrostratus (Cs), 
and nimbostratus (Ns). With this scheme, more than one of the observed and 
recorded cloud types may fall into one of the nine basic cloud categories as 
delineated above. This is especially true in the ship reports where a total of 
27 cloud types are identified. The conversion of the cloud types on tape to the 
nine basic cloud categories was achieved in the processing by incorporating the 
conversion tables, shown in Tables 4-1 and 4-2, into the programs. In examin- 
ing the data contained in these tables, a number of objections must be dealt with. 

TABLE 4-1 

CLOUD-TYPE CATEGORIES FOR LAND STATIONS 


Tape Code 

Defi nition 

Cloud Type 

0 

None 


1 

Fog 


2 

Stratus 

St 

3 

Stratocumulus 

Sc 

4 

Cumulus 

Cu 

5 

Cumulonimbus 

Cb 

6 

Altostratus 

As 

7 

Altocumulus 

A c 

8 

Cirrus 

Ci 

9 

Cirrostratus 

Cs 

K 

Stratus Fractus/Fractostratus 

St 

M 

Cumulus Fractus/practocumulus 

Sc 

N 

Cumulonimbus Mamma 

Cb 

O 

Nimbostratus 

Ns 

P 

Altocumulus Castellanus 

Ac 

R 

Obscuring Phenomenon 


r t _ n 

Obscuring Phenomenon other than Fog 
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The most serious, perhaps, is that "fog," which in essence is stratus extending to 
the surface, is not accounted for in the cloud statistics. This decision was made 
necessary because the unconditional cloud cover statistics, previously derived by 
Sherr et al, do not account for the presence of "fog" as legitimate cloud cover. 

Since one of the primary goals is to establish a consistent set of cloud data, it was 
decided not to incorporate "fog" into the cloud -type statistics. For most regions, 
the exclusion of "fog" is not expected to introduce serious errors in simulation. 

The regions most likely to be affected are in the northern latitudes. 

A second deficiency in the compilation of cloud types is that in processing 
the ship data the cloud category Nimbostratus (Ns) was not used. This was made 
necessary because Nimbostratus was not used as an explicit cloud group by the 
ship reports. Apparently, when it did appear it was reported as a multilayered 
cloud. If the statistics from the other stations are any indication, however, it 
does not appear that Nimbostratus, as a separate cloud group, occurs with suf- 
ficient frequency to cause any real problems in simulation. Other objections to 
the cloud groups selected for processing undoubtedly exist. However, it should 
be emphasized that the primary purpose of the statistics is not for detailed clima- 
tological studies but for the practical problem of cloud attenuation simulations for 
earth observations. For this latter purpose, no detailed breakdown of cloud types 
is really essential since the microphysical properties of many of the cloud types 
are not essentially different. 

The problem of compiling cloud-layer statistics was slightly simplified by 
assuming that no more than three cloud layers can occur at any given observation. 
While it is true that situations with four or more cloud layers are sometimes 
reported, the frequency of occurrence is so low that to provide stable statistical 
samples would require extremely long record periods. Instead, the computer 
programs were designed to treat all occasions of three or more cloud layers as 
three cloud layers. Statistics were therefore compiled only for 1, 2, and 3 cloud 
layers . 

It was previously noted in Section 3 that for a number of stations, cloud-type 
data are recorded on tape for the standard GMT observation times. The computer 
programs were designed to substitute the required local time closest to the observed 
Greenwich Mean Time observations. Basically, this procedure implicitly assumes 
that no significant changes in cloud-type distribution have occurred during the hours 
between the true observation time and the observation time now associated with 
the cloud statistics (see Table 3-6). The validity of this assumption is obviously 
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dependent on the cloud type and the region. For a cloud type which exhibits marked 
diurnal dependence, this extrapolation may not be very good for certain periods of 
the day. However, it should again be emphasized that the primary purpose of the 
statistics is to provide practical data for simulation purposes and not a definitive 
global cloud climatology. 

In addition to the above statistical parameters, the programs were also able 
to compile the unconditional cloud cover statistics in the five basic categories used 
by Sherr et al and defined previously in Table 2-1. A number of safeguards were 
built into the program to recognize and reject observations containing obscurations 
and wrong codes, and to handle situations with missing data. The products of the 
processing were printouts of the compiled statistics in two formats and punched 
cards on which the statistical information is stored for future use in simulation 
programs. The formats of these outputs are discussed in Section 4.3. 

4.2.2 Microfilm Data 

The nine basic cloud types and three cloud layer groups used in processing 
the taped data were also employed in compiling statistics from the microfilm data. 
However, because of the time and effort involved in extracting the data by hand, no 
attempt was made to derive individual sets of cloud-type and cloud-layer statistics 
for each of the five cloud cover categories. Instead the assumption was made that 
for these stations and regions, the statistics are independent of cloud cover so that 
only a single set of cloud-type and cloud-layer statistics need be derived instead of 
four. In addition, the unconditional cloud cover statistics were not explicitly derived. 
However, with the data, statistics for the clear-sky category can be derived from the 
difference between the sum of the number of cloud layer observations, a value obtain- 
able from the cloud layer statistics, and the total expected observations, a value 
derivable from the length of the record. In fact, this procedure was adopted in 
deriving the statistics for cloud cover category 1 (clear skies) compiled from the 
microfilm records to compare with the corresponding statistics derived by Sherr 
et al. 

The problem of time extrapolation encountered in processing the taped data 
was much more acute when it came to reducing the microfilm data. In the process- 
ing of the taped data, the difference between actual observation time and the assiimed 
time recorded on punched cards never exceeded an hour. With some of the micro- 
film data, this time difference was sometimes as large as four hours. These dis- 
crepancies arise because, for a few of the stations, no observations are available 
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for the early morning and late evening hours. Observations at Pelly Bay are an 
example. The earliest daily observations at this station were made at 05 LST and 
the last observations were at 17 LST. Since no other data are readily available, 
the observations made at these times were simply extrapolated to 01 LST and 13 
LST. Similar extrapolations were made with the other off— time observations. 

4.3 Data Formats 

4.3.1 Taped Data 

Since the statistics are to be used primarily in computer simulations, the 
data were preserved on punched cards. In past studies, the practice has been to 
store on cards the actual derived ’’probability*’ values. This practice had the 
virtue of providing the cloud statistics in a form directly usable in simulation 
procedures. On the other hand, it has been found that by storing only the norma- 
lized values, it is impossible to improve on the statistics through the simple pro- 
cedure of enlarging the data sample since information on the size of the original 
data sample is not preserved. This latter consideration was taken into account in 
the design of the punched card format shown in Figure 4-1, In this format, two 
cards will be used to code the data for the statistics from a single hour. With 
eight three-hourly sets of statistics per day, one month’ s statistical data from a 
given region can be punched onto 16 cards. In recording the data, the following 
will be used: 

Columns Comments 

1-2 2 digit month identification code with 01 = Jan 12 December. 

3-4 2 digit cloud cover region identification code 

5 1 digit local time (LST) identification code with 1 = 0100 LST, .... 

8 = 2200 LST 

6 1 digit cloud amount (”j”) code with the following values; 

For ”j” = 2; all data in columns 7-42 apply to cloud cover 

category 2, and all data in columns 43-72 apply 
to cloud cover category 3. 

For ”j" = 4; all data in columns 43-72 apply to cloud cover 
category 5. 

7 - 33 3 digit (000 to 999) frequency count, f(j,i), of the occurrences of each 

of the 9 cloud groups for cloud cover category 2 or 4. 

34 - 42 3 digit (000 to 999) frequency count, f(j,k), of each of the 3 cloud layer 

groups for cloud cover category 2 or 4. 
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Figure 4-1 Punch Card Format 


Columns 


Comments 


43 - 69 Same as columns 7-33 but for cloud cover category 3 or 5. 

70 - 7S Same as coiuiims 34-42 but for cloud cover category 3 or 3. 

79 - 80 The four column locations of the two-card set for each of the three 

hourly sequence will be used together to enter the three-digit fre- 
quency count, the event of "cloud cover category 1,’* i.e., f(j=l). 


In order to provide usable data for publication, printouts formatted to 
display the data were obtained. Figure 4-2 shows an example of the printout of 
one month* s data (Month 5) for one region (Region 1). The cloud-type and layer 
data for each of the indicated Local Standard Time (LST) periods are contained 
in two matrices. The first matrix, with 4x9 elements, contains the frequency 
counts of each of the cloud types indicated (along a column) for each of the cloud 
cover categories shown (along a row). The second matrix, with 5x3 elements, 
tabulates the cloud layer frequencies for the five cloud cover categories and the 
three cloud layers shown. The last number, to the right of the cloud layer stati- 
stics matrix,^ is the frequency of the clear skies (or cloud cover category 1). 

A printout in the form shown in Figure 4-2 was obtained for each month 
and for each region. These are reproduced in Appendix A in the reduced size 
format shown in Figure 4-3. With these data, it is possible to reproduce the 
card deck for further processing and analysis. 

In both data formats, the punched cards and the printouts in Appendix B, 
only the actual sample frequencies are preserved and not the probabilities. How- 
ever, it can be shown that all necessary probabilities can be readily obtained 
from the tabulated data. Consider, for example, the probability P(j,i) of the 
occurrence of cloud type i when it is known that the cloud cover category j has 
occurred. This can be computed simply from the tabulation by: 


where 


P(j»i) 


i=9 

S f(j.i) 
i = l 


(4-1) 


f(j,i) is the tabulated frequency of occurrence of cloud type i 
when the cloud cover category is j. 


Similarly, the probability, P(j,k), of the occurrence of k cloud layers, when the 
cloud cover category is known to be j, is given by 


P(j.k) = • (4-2) 

S f(j.k) 
k=l 
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Figure 4-3 Format for Appendix B 
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If, on the other hand, the probability P(i or k) of occurrence of cloud type i, or 
k cloud layers, is required independent of the cloud cover j, it can be computed 
from 


P{i or k) = 


f(j = TOT, i or k) 


2 

i or k 


f(j = TOT, i or k) 


where 


5 

f(j = TOT, i or k) = r f(j, i or k), 

j=2 


(4-3) 


(4-4) 


Furthermore, with the frequency count of category 1 for the unconditional cloud 
cover statistics being provided on the punched cards and printouts, it is even 
possible to compute the probabilities P(j) of all the unconditional cloud cover 
categories by 


P(j) = ■■ ■ (4-5) 

S f(j) 

3 = 1 


where 


from 


f(j) is the frequency of occurrence of cloud cover category j computed 


3 

f(j) = Z f(j,k) for j 1. (4-6) 

k=l 

f(j,k) is, as before, the frequency of occurrence of k cloud layers when the cloud 
cover category is j. 

Once programmed to perform the required summations and divisions, the 
computer processing time to generate the required probabilities using the equali- 
ties identified is minimal. In providing the statistics for the analysis discussed 
in the next section, the data were in fact processed in such a manner. Figure 
4-4 shows an example of the data format used. For each joint event, both the 
frequency count and the probability were displayed. 


31 



32 
































4.3.2 Microfilm Data 


The statistics derived from the microfilm data are recorded on cards and 
in Appendix A in the formats previously described for recording and displaying 
the data derived from the tapes. The cloud-type and cloud-layer statistics are, 
as before, separately displayed and tabulated for each of the four nonclear-sky 
cloud cover categories. However, close examination of the data will show that 
these statistics do not change with cloud cover category, since in fact only one 
set of these statistics were derived. The repetition of these statistics through 
the four cloud cover categories was made simply to provide a uniform data 
format for all stations. In the context of the discussions in Section 4.3. 1, the 
statistics tabulated are in fact the values of 

5 

f (j = TOT, i or k) = Z f(j or k). 

j=2 

A second difference in the data tabulations for the microfilm data is the 
lack of statistics for the unconditional clear-sky frequencies. While these data 
can be derived from the tabulaietl layei frequcncicc as previously rioted 

they are of no real value in deriving the statistics of the remaining cloud cover 
categories by Eqs. (4-5) and (4-6). For this reason, the statistics were not 
included. The only available sets of unconditional cloud cover statistics for these 
regions are those developed by Sherr et al. 
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5. DISCUSSION OF THE DERIVED CLOUD STATISTICS 


The compilation of cloud type and layer statistics has produced a wealth of 
rxteresting climatological results. In this section, some of these are discussed. 
However, the discussion will be restricted to only those features of the statistics 
Df significance to the problem of cloud simulation for earth observations. In parti- 
cular the questions of homogeneity of regions, sample size, diurnal and annual 
/ariabilities, and regional differences are examined. 

5. 1 Regional Homogeneity 

In Sherr et al (1968), the question of regional homogeneity was examined 
In detail with regard to the unconditional cloud cover distributions. It was shown 
then that the differences in the unconditional cloud cover statistics between stations 
in the same region were insignificant for some of the regions. Of the six stations, 
from Minot, North Dakota, to Furth, Germany, examined for Region 11, for example, 
the differences in the percentage frequency of occurrence for any given cloud cover 
category were generally less than 5%. However, other regions showed less homo- 
geneity. A case in point is Region 3, Foi the five stations analyzed for this 
region, the differences in the percentage frequency of occurrence for the cloud 
cover category 5 (total sky cover) in some instances exceeded 40%. While it 
is not the intention of this study to improve on the homogeneity of the regions 
by redefining regional boundaries, the limitations of the data introduced by 
what are, at best, quasi-homogeneous regions should be discussed and under- 
stood. By so doing, it is hoped that the data developed in this study and in the 
previous study will not be misused in applications for which they were not 
intended. 

The question of regional homogeneity has increased significance when 
the cloud type and layer statistics, developed in the present study, are used 
jointly with the unconditional cloud cover statistics previously developed by 
Sherr et al. While the cloud climatological regions remained unchanged, the data 
samples from which the two sets of statistics were derived are quite different. In 
most instances, this difference resides in the length of data record used for statis- 
tical compilations. However, more serious are the changes in '’prototype” stations 
for some of the regions in the current study, made necessary by the reasons cited 
previously in Section 3. If the cloud climatological regimes are quite different for 
the two stations from which the two sets of cloud statistics were derived for the 
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same region, the statistics would not be, in fact, self consistent, thus diminish- 
ing the quality of the entire data set. A study was therefore conducted to deter- 
mine the similarity in cloud climatologies between two stations, both used as 
prototype stations, for the same region. 

It would be highly desirable to examine the homogeneity of the cloud regions 
with respect to cloud type as well as cloud cover. However, with the data available 
to this study, this was not possible except for Region 12. For this region, micro- 
film WBAN data of cloud cover and cloud type were available and analyzed for Bien 
Hoa and Ton Son Nhut, both stations located near Saigon. Because these stations 
are located no more than 3 0 km apart, the comparison between their corresponding 
cloud-type statistics provides a measure of the ’’best” one should expect in cloud 
region homogeneity. As for the other regions, the only type of statistics available 
for comparison are the distributions of unconditional cloud cover categories. It is 
here assumed that the degree of similarity in the cloud cover statistics between 
two stations is a measure of the degree of similarity in the cloud-type distributions 
between the stations. 

From Table 3-4, it may be seen that changes in the prototype or represen- 
tative stations were made for Regions 9, 10, 12, 13 ond 14. Comparisons cf the 
unconditional cloud cover statistics derived previously by Sherr et al with those 
derived in this study for Regions 9, 13 and 14 are shown in Figures 5-1 to 5-3. 

These figures, following the format used in Sherr et al, illustrate the probability or 
frequency of occurrence of each of the five cloud cover categories for the months 
and times indicated. Similar figures cannot be prepared for the other two regions, 
Regions 10 and 12, since the cloud cover statistics by categories were not extracted 
for the new representative stations. Felly Bay and Ton Son Nhut. Instead, the 
annual variability of the important category 1 (clear sky) statistics was derived from 
the cloud layer statistics for each of these two stations. The results are shown 
plotted in Figures 5-4 and 5-5. Figures 5-6 and 5-7 have been included to illus - 
trate similar annual variabilities for Regions 13 and 14. In the latter two figures, 
statistics for cloud cover category 5 (total sky cover) are also shown. 

The data presented in this sequence of figures show that the differences in 
cloud cover statistics between the pairs of stations are well within the limits pre- 
viously found by Sherr et al. Even in those instances where large departures are 
apparent, as in the case of Region 10, the large differences are found only during 
one or two months of the year. In this latter case, the surprising fact is that 
despite the great differences in geographical locations between Harbin, China, and 
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Figure 5-2(a) Comparison of Region 13 Cloud Cover Statistics 
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Figure 5-3{b) Region 14 Cloud Cover Statistics 
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Figure 5-4 Region 10 Cloud Cover Statistics 
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Figure 5-5 Region 12 Cloud Cover Statistics 
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Figure 5-6 Region 13 Cloud Cover Statistics 
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Figure 5-7 Region 14 Cloud Cover Statistics 
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Pelly Bay in the N.W.T., the shapes of the curve defining annual variations of 
dear-sky probability, with their characteristic minima during the summer 
months, are in fact very similar. 

The microfilm data from Vietnam afforded an opportunity to examine the 
differences in the cloud statistics between two stations which are located close to 
each other. These tv*o stations, Bien Hoa and Ton Son Nhut, are approximately 
30 km apart. In comparing the data from these two stations, it was found that the 
distributions of cloud-type statistics were on the whole very similar as expected. 
However, even for stations so close to each other, some instances of significant 
difference in the frequency of occurrence for a particular cloud type were found, 
as may be seen in Figures 5-8 and 5-9. Particular differences of note are the 
statistics for cloud type 2 (stratus) 01 LST - January (Figure 5-8); the overall 
distributions for 01 LST - April (Figure 5-8); and the distributions of cloud types 
2 (stratus) and 7 (cirrus) for 01 LST - July (Figure 5-9 ), While it is not intended 
that these data be interpreted as being representative of the expected in-region 
differences in cloud-type statistics for all regions, they do illustrate that a high 
degree of homogeneity in statistics for any region should not be expected. 

5.2 Sample Size Dependence 

In previous studies, no analysis was made of the dependence of cloud stati- 
stics on size of the sample. The implicit assumption was that the larger the size 
of the sample, the more representative are the derived statistics. This is certainly 
a valid assumption and with the long periods of records available in the form of 
cloud-cover summaries, the influence of the sample size on the statistics was never 
a real question. In this study, a limited opportunity was afforded to examine the 
dependence on sample size of not only cloud cover statistics but also of the cloud 
type and cloud layer statistics. This was made possible by the data-proces sing 
adopted in the initial phases of this study. In preliminary tests of the programs 
to extract relevant statistical data from the taped WBAN records, for example, 
printouts of the statistical data were obtained for a few of the stations for record 
lengths shorter than the maximum available on the tapes. Furthermore, for 
those stations for which only microfilm data were available, the statistics for 
individual years were separately tabulated before being combined together to form 
the final data sample. Some of these results are presented here purely as examples 
to illustrate possible influence of sample size on the statistical distributions. 

An example of the comparison of the statistics derived from two sample 
sizes is shown in Figure 5-10. Here the unconditional cloud cover statistics 
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Figure 5-9 Region 12 Statistics 
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Figure 5-10 (a) Cloud Cover Sample Size Comparison 
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Figure 5-10(b) Cloud Cover Sample Size Comparison 
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derived from 5 years and 10 years of taped recorded observations for Region 16 
and for the months and times shown are plotted. While the match between the two 
sets of corresponding data is not perfect, it is clear that for all practical purposes, 
the differences are negligible. 

A more complete set of data were derived from Region 5 (San Diego). The 
taped data for this region contained a full 10 years » record of WEAK observations. 
These data were processed to produce statistics based on 1-year, 5-year, and 10- 
year records. Figures 5-11 to 5-13 present representative comparisons for the 
three sample sizes. The data were selected to show these comparisons for January 
and July, and for 01 LST and 13 LST . The first set of data, in Figure 5-11, shows 
the comparisons for unconditional cloud cover statistics. The second set. Figure 
5-12, shows the dependence of cloud-type statistics on sample size. The third set, 
Figure 5-13, shows the data for cloud layer statistics. In examining these figures, 
it IS clear that at least for this region, the unconditional cloud cover statistics are 
the least sensitive to sample size. Furthermore, for these statistics, the differ- 
ence between the data derived using the 5 -year sample and the 10-year sample is 
for all intents and purposes negligible. 

The cloud-type statistics show more dependence on sample size. This is 
particularly evident in comparing the data derived from the 1-year sample with 
those derived from the larger sample size. The ’’amplitudes” in the 1-year data 
appear exaggerated such that the dominant cloud types have a too large probability 
of occurrence while the probabilities of occurrence for the remaining cloud types 
are much too low. This characteristic of the small sample statistics is undoubtedly 
due to the distribution characteristic of the particular year selected for process- 
ing. However, it is clear from the figure that the statistics derived from sample 
sizes of 5 years or more records approach stable values. 

The greatest variability is in the cloud layer statistics. While, in general, 
the statistics favor cloud cover situations consisting of a single layer, there appears 
to be a tendency to increase the probability of occurence of cloud cover situations 
with multiple cloud layers with increasing sample size. In comparing all of the 
data, it does appear that overall the differences between the 5-year and 10-year 
data are much smaller than the differences between the 1-year sample and data 
derived from any of the two larger samples. However, exceptions are in evidence, 
as shown in the data for Month 7, 13 LST (Fig. 5-1 3c) and for cloud cover category 
4. In this case, the trend in the data derived from the 5 -year data sample is total- 
ly unlike that in the data derived from the 1 -year and from the 10-year data samples. 

In order to examine further the dependence of cloud-type statistics on 
sample size, the data from Region 12 (Ton Son Nhut) were analyzed in detail. 
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Figure 5-1 1(a) Region 5 Cloud Cover Sannple Size Comparison 
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Figure 5-13(a) Region 5 Cloud Layer Sample Size Comparison 
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Figure 5-1 3(c) Region 5 Cloud Layer Sample Size Comparison 
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Data from Ton Son Nhut were readily available for the detailed analysis, and the 
climatology of this region is such that the frequency of occurrence of cloudy skies 
is quite high. Under this condition, there is a likelihood that the cloud-type 
statistics would be quite stable from sample to sample even for sample sizes 
as small as those derived on the basis of a single year of observation. This a 
priori conclusion appears to be well borne out by the data shown in Figure 5-14. 
Here, for the three years of 1967, 1966 and 1969, the diurnal changes in the oc- 
currence of some of the major cloud types, Cu, Cb, Ac, and Ci, are shown plotted 
for the months of January and July. As shown in the figure, the year-to-year 
differences in the cloud -type frequency distributions are small for this region as 
anticipated. 

It might be recalled that the unconditional cloud cover statistics derived by 
Sherr et al were based on data obtained over long periods of, records, generally 
in excess of 15 years. For Regions 3 and 14, the statistics were derived from 
Revised Uniform Sumnnaries of observations obtained over periods of 21 years 
and 25 years, respectively. Samples of the unconditional cloud cover statistics 
based these records are shown in Figures 5-15 and 5-16. Also shown are the 
corresponding statistics derived from the smaller data sample (6 years for Region 
3 and 10 years for Region 14). Again, the improvement” in the statistical distri- 
butions resulting from the longer data samples are not significant. 

5,3 Regional Variability 

Using the cloud -type data generated, a study was conducted of the differ- 
ences in these statistics between the regions. This study was limited to a com- 
parison of statistics between pairs of regions in the low latitudes, the midlati- 
tudes, and in the high latitudes, 

5. 3. 1 Low Latitude Figures Regions 

The pair of low latitude regions selected for comparison was Region 1 
(Dhahran, Saudi Arabia), and Region 3 (Angeles, Luzon, P. I,), According to 
the map of the cloud regions shown in Figure 2-1, Region 1 is located in major 
desert areas while Region 3 is in more of a tropical oceanic regime. Sherr et al 
prescribed no major cloud type for Region 1 but did indicate convective type clouds 
to be the predominant type of clouds for Region 3, based purely on climatological 
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Figure 5-14(a) Region 12 Cloud Type Frequency Distributions 
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Figure 5-14(b) Region 12 Cloud Type Frequency Distributions 
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Figure 5-15 Region 3 Sample Size Comparisons 
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Figure 5-16 Region 14 Sample Size Comparisons 
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considerations. Examples of the actual cloud-type distributions for these two 
tropical regions are shown in Figures 5-17a and 5-17b, and 5-18a and 5-18b. As 
expected) these figures show that despite the latitudinal similarity in the locations 
of the regions# the distribution of cloud types is quite different. Also, as expected 
from the maritime climate of Region 3, the cloud types are predominantly convec- 
tive with cumulus in the afternoon (13 LST) and stratocumulus in the very early 
morning hours (01 LST). More importantly, these clouds are formed locally in 
the presence of a sufficient supply of water vapor. The predominant cloud types 
for Region 1, on the other hand, are middle and high clouds. These clouds are 
most likely the result of spreading from synoptic or even planetary scale events 
rather than any local effect. This is especially true in the case of cirrus which 
appears to have an overall higher probability of occurrence than any of the other 
cloud types. The differences in the cloud type distributions between these two 
stations clearly demonstrates the influence of local effects in the climatology of 
clouds . 


5. 3. 2 Nlidlatitude Regions 

In the midlatitudes. Regions 4 and 8 were selected for comparison. Region 
4, with its representative station, Tampa, Florida, is still basically subtropical 
in nature, subjected to the local formation of convective clouds. Region 8, with 
Moimtain Home, Idaho, as representative station, is clearly midlatitude. Accord- 
ing to Sherr et al, this latter region should have cloud types primarily associated 
with synoptic scale systems. Based on these simple considerations, the cloud 
type distributions for the two regions can be expected to be quite different. 

Figures 5-19a, 5-19b, 5-20a and 5-20bshow examples of the cloud-type 
probabilities for these regions. These data in fact show that despite the differences 
in the general climatic regimes, the cloud-type probability distributions for the 
two regions are remarkably similar. The major difference, as expected, lies in 
the more probable or frequent occurrence of cumulus for Region 4. However, the 
magnitude of the difference is less than expected. 

The reason for this apparent similarity in the cloud -type statistics is 
perhaps because Region 4, while subtropical in nature, lies sufficiently close to 
the tracks of cyclonic storm systems to have its cloud statistics influenced by 
their passage. Hence, the frequent occurrence of middle layered clouds. 
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Figure 5-17(a) Low latitude Region Comparison 
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Figure 5 -18(a) Low Latitude Region Comparison 
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5. 3. 3 High JLatitude Regions 

For the high latitude comparison, the data from Region 10 (Pelly Bay) and 
Region 15 (Resolute) were used. Examples of these data are shown in Figure 5-21. 
Both of these regions were identified by Sherr et al as having predominantly clouds 
associated with synoptic scale systems. The data actually shows that low strati- 
form clouds appear to be the dominant cloud type. However, these statistics may 
be deceptive. It is well known that low stratiform clouds often obscure the pres- 
ence of upper cloud layers. Therefore, even though the data may show low clouds 
to be frequently present, it does not necessarily imply that middle and high clouds 
do not occur as often. This problem of high and middle cloud layer obscuration by 
low clouds is not just confined to these high -latitude stations. Statistics for the 
high and middle cloud types for any region can be expected to be on the low side 
when the basic data used are derived from ground observations. For most regions, 
the difference between the true statistics and those derived from ground observa- 
tions are not expected to be of great significance. For regions with frequent oc- 
currences of stratus and cloud cover category 5, however, this problem may be 
acute. 
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6. STATISTICS OF IN-CLOUD PARAMETERS 


With the global statistics of cloud type and cloud layer discussed in Section 
5, the establishment of global statistics of the four -dimensional distributions of 
cloud drops and liquid water contents requires only models of the intensive para- 
meters appropriate to each of the nine cloud types. The inclusion of cloud thick- 
ness and height of cloud base, would permit a vertical profiling of the microphysi- 
cal cloud parameters. In this section, the research to establish these model 
parameters and the model parameters themselves are discussed. 


6. 1 Analytical Considerations 

Measurements of in-cloud parameters are at best difficult to obtain as 
attested by Mason (1962) among others. To make maximum use of the data which 
are available, empirical equations have been derived to analytically fit the meas- 
urements and to provide a means to interrelate the voiious parameters. These 
equations reduce the modeling of the microphysical structures of clouds to the 
specifications of a small set of parameters. In the present context of providing 
a cloud data bank ior computer simulation this property of the analytic approach 
has great merit and was investigated. 

Within any volume in a cloud, there is a spectrum of droplet sizes. 
Fortunately, theoretical considerations suggest, and measurements prove, that 
the distribution of droplets among the various droplet sizes is systematic, and 
therefore amenable to analytic expressions. Measurements by Diem (1948), 

Squire (1958) and others all show that droplet size distributions are generally 
bell-shaped, and more or less asymmetrical with a tail extending towards larger 
radii. It is readily seen that such distributions can best be fitted by an exponential 
function. A number of empirical formulas, differing only in minor details, have 
in fact been proposed, and in many instances used to specify cloud droplet spectre . 
In the research conducted during this study, the formulation discussed by A. M. 
Borovikov et al (1961) is perhaps most complete. 

The expression for the size spectrxim of cloud droplets is given by 

n(r) = ar^ exp (-br) (6-1) 

where 

n(r) is the volume concentration at radius r (#/cm - p) 

r is the droplet radius 
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and 


<x <i o are constants adjustable to fit the particular spectrum 
in question. 

The droplet volume density or number concentration N (r?/cm^) is then given by 


oo 


N 


n(r )dr . 


( 6 - 2 ) 


Using Eq. (6-1) for n(r), the integral can be equated in closed form as 


N 


J [ a r"" exp {- br)] dr - 
« b 


(6-3) 


Furthermore, it can be shown that the mean radius, r , is given by 


00 


* r n( 


r) dr 


oo 

J n(r) dr 


3 

b 


(6-4) 


And by differentiating Eq , (6-1) with respect to r, and setting the result equal 
to zero, the modal radius, r^, of the distribution is obtained as 

^d ^ f- (6-5) 


Since the mean cloud liquid water content W 


igm/m) is 


expressed as 


W 


4 

3 


(10) ^ 7T r^ N, 


( 6 - 6 ) 


W may also be expressed in terms of the constants a and b, by substituting 

Eqs. (6-2) and (6-4) for N and r , as 

m 

w = 10'^ a , 

b^ 

It is clear from the above discussion that by using Eq. (6-1) as the empirical 
equation for drop size distribution, an internally consistent model of in-cloud 
parameters is achieved by specifications of just r^ (or r^) and W, and the geo- 
metric parameters of cloud thickness (T) and cloud base altitude (H). 
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It should be pointed out that the above formulation assumes that the droplet 
size spectrum can be fitted by Eq. (•!)• Examination of actual spectre will show 
that no spectrum is really continuous or smooth. Under these conditions, no 
analytic expression can be expected to fit perfectly all aspects of even a single 
droplet spectrum much less to fit all spectra. For the purposes of determining 
approximate evaluations of order -of -magnitude cloud -attenuating effects for en- 
gineering purposes, the difference between the actual values of the in -cloud para- 
meter and those determined on the basis of the analytic techniques are expected to 
be inisgnificant. However, if computations of the detailed optical properties of 
clouds are required, the analytic expressions are perhaps inadequate. 


6.2 Data Sources 


In-cloud and geometric parameters of clouds are not routinely observed 
and reported as are cloud types. The only possible exception is cloud base 
altitude (H), This parameter is observed ....d reported part of the surface 
observation procedures performed at weather observation stations. However, 
the height of cloud bases, especially for high and middle clouds, are often esti- 
mated (eyeballed) without the benefit of objective instrumentation. Consequently, 
there is no large body of empirical observations which could serve as the raw 
data base from which meaningful statistics of the global distributions of r^ (or r^). 
W. T. and H. can be derived for each of the cloud types. The only data available, 
are sets of experimental measurements made of a limited number of locations, 
and at a limited number of altitudes. To further compound the problem, most of 
the different sets of experimental measurements were made generally using dif- 
ferent types of sensors, employing both direct and indirect measurement tech - 
niques. In some cases, the equipment itself was experimental. Indeed, even now. 
neither the measurement techniques nor equipment are standardized. It is there- 
fore not possible to provide a definitive set of in-cloud parameters at this time. 

In reviewing the literature on the microphysical structure of clouds, it 
becomes abundantly clear that while the Uterature itself is quite extensive, the 
actual data upon which the literature is based is extremely limited. Repeated 
references are made to the now classic measurements made by Diem (1948). 
Neuburger (1949), and Weickman and Aufm Kampe (1953) among others. 
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Other papers either deal with the application and generalization of these measure- 
ments, or reinterpret these measurements in summary form. This latter cate- 
gory of papers and publications provide the best source of data for the establish- 
ment of cloud model parameter In the literature review, it was ascertained 
that the best of these summaries is the work on cloud physics by Borovikov 
c‘t al (I96lj, This publication not only summarizes the quite extensive measure- 
ments made in the Soviet Union, but the data obtained by Western European and 
Ame rican worker s are also well represented. More importantly, this publication 
discusses in detail the variability of in-cloud and geometric parameters using all 
data available. However, for the purpose of providing a representative tabulation 
of significant cloud parameters, the summary by Gaut and Reifenstein (1971) in 
their study of microwave radiation interactions with the atmosphere was found to 
be most appropriate. The tabulated data are based on the works of a number of 
investigators in the United States and Great Britain. These include: Valley (1965), 
Fletcher (1966), Mason (1957), Blau et al (1966), Deirmendjian (1964), and Crane 
(1966), 

6. 3 Tabulation of Parameters and Discussion 

In the study by Gaut and Reifenstein previously referenced, a total of nine- 
teen cloud and precipitation models are defined. Each of these models is specified 
by a set of parameters including base and top altitudes, density or liquid water 
content and two shape parameters. Except for the latter parameters, made neces- 
sary by Dei rmendjian’ s (1964) approach to the analytic representation of the drop 
size spectrum which they adopted, their model parameters, with some minor 
modifications, are directly applicable to the problem in hand. Table 6-1, the 
tabulation of significant in-cloud and geometric properties, was based on a simpli- 
fication of the cloud models defined by Gaut and Reifenstein (1971). 

The values shown in Table 6-1 are meant only to be representative and not 
definitive. Indeed, the literature suggests that any attempt to produce definitive 
cloud models would be futile. As pointed out in Borovikov (1961), the cloud is 
never stationary, and in fact "... Droplets form, grow, and then either fall out 
of the cloud or again shrink by evaporation. In the first stage of development the 
cloud is composed of fine droplets and is relatively monodispe rse ; subsequently 
the droplets grow larger and their spectrum expands. Finally, in the last stage, 
the "dissolution" of the cloud may take place either as a result of evaporation of 
its droplets or following precipitation and later evaporation of the remainder of 
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TABLE 6-1 


MODELS OF CLOUD PARAMETERS 


Cloud Type 

Base 

(m) 

Top 

(m) 

Modal 
Radius 
(rm, (p) 

Density or Liquid 
Water Content 
{g/m}) 


500 

1, 000 

10. 0 

0.50 

Cumulus 

1, 000 

1. 500 

10. 0 

1.00 


1.500 

2, 000 

10. 0 

0.50 

St r at o c umulu s 

500 

1, 000 

10. 0 

0.25 

Stratus 

150 

1, 000 

10. 0 

0. 25 


0 

300 

400. 0 

6.30 


300 

1, OUU 

20. 0 

7. 00 


1, 000 

4, 000 

10. 0 

8. 00 

Cumulonimbus 

4, 000 

6, 000 

10. 0 

4. 00 


6, 000 

6, 000 

10. 0 

3. 00 


8, 000 

10, 000 

40. 0 

0.20 

Altostratus 

2,400 

2, 900 

10. 0 

0. 15 

Altocvimulus 

2,400 

2.900 

10. 0 

0. 15 

Cirrus - Arctic 

WBM 


40. 0 

0. 10 

- Mid-JLat, 

WKSM 


40. 0 

0. 10 

- Tropical 

IB 


40. 0 

0. 10 

Cirrostratus - Arctic 


6. 000 


0. 10 

- Mid-Lat* 


7, 000 

BB9 

0. 10 

- Tropical 


8, 000 

40. 0 

0. 10 


0 

500 

200. 0 

1.00 

Nimbo stratus 

500 

1,000 

10. 0 

2. 00 


1, 000 

2,000 

10. 0 

3.00 


2,000 

4, 000 

10.0 

2.00 
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the cloud," As the spectrum, n{r), and total number density, N, change, the 
liquid water content, W, also changes. Furthermore, the detailed development 
of the droplet spectrum is a function of the height level within the cloud. Figure 
6-1 shown examples of this height dependence. The data were taken from the work 
by Borovikov (1961) and show the dependence of the mean droplet radius, r , on 
the location of the sampled volume within stratiform clouds. It is clear from this 
figure that apart from the vertical variability in the greater variability ap- 

pears to be between samples. Measurements of liquid water content also exhibit 
larger variabilities as clearly exhibited in Figures 6-2 and 6-3. The data are 
again taken from Borovikov. Because of these empirical results, the in-cloud 
parameters of a few of the models shown in Table 6-1 are separately tabulated 
for a number of layers extending from base to cloud top. 

Measurements suggest that the geometric parameters of cloud thickness, 

T, and cloud base altitude, H, can be highly dependent on latitude. Figure 6-4, 
based on measurements compiled in Borovikov, shows the strong latitudinal 
dependence of cloud base altitudes for stratocumulus and a very weak dependence 
for stratus. However, there is really insufficient data to define the degree of 
dependence. For this reason, the data in Table 6-1 show latitudinal dependence 
only for the two high cloud types. The dependence is based purely on the assump- 
tion that cirrus and cirrostratus are formed near the tropopause which, itself, 
decreases in altitude from the Equator to the Poles. 

6.4 Use of In-Cloud Parameters in Simulation 

It is at present proposed that the in-cloud and geometric properties modeled 
in Table 6-1 be used in simulation of the vertical distribution of cloud droplets. 
Obviously, as presented, the variabilities of these parameters are not incorporated 
into the models, and therefore not simulated. One possible approach to provide 
practical data to incorporate in-cloud variabilities in simulation is to assume 
that the statistical distributions of the cloud parameters are normal or Gaussian. 
The cumulative probabilities for each of the parameters can then be computed 
from standard probability tables, and the results used in a Monte Carlo scheme 
for the selection of the appropriate set of in-cloud parameters. However, 
in order to fully implement this approach, it will be necessary to include in the 
cloud models the standard deviation or variance of each of the parameters. These 
values, together with the mean value will uniquely specify the Gaussian distribu- 
tion corresponding to each of the parameters. 
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Ctoud Type - Stratocumulus 



Height b Cloud (H= Total Cloud Thickness) 

Figure 6-1 (a) Stratocumulus Droplet Variation 
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Cloud Type • Stratus 



Height in Cloud (H= Total Cloud Thickness) 


Figure 6-l(b) Stratus Droplet Variation 


82 





Cbud Types > Stratocumu/us 
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Figure 6-3 Cloud Water Content 


84 


E 



0 ' 10 20 30 40 50 60 70 

LATITUDE DEGREES 


I I I 



0 10 20 30 40 50 60 70 


LATITUDE (degrees) 

Figure 6-4 Fatitudinai Dependence of Cloud Base Altitudes 
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This suggested apprud.ch wiii add sophistication to the modeling technique. 
However, in view of the total lack of systematic data to adequately establish any 
idea of the population variances, we do not feel that any real improvement in the 
simulation will result from its implementation. 
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7. SIMULATION: A SUGGESTED APPROACH 


In Section 4, a possible approach to a simulation scheme using cloud-type 
statistics was briefly suggested. In this section, this approach is discussed in 
greater detail. The suggested scheme is primarily an extension of the Monte Carlo 
technique originally proposed by Sherr et al for the cloud simulation problem using 
cloud cover statistics. For this reason, a summary recapitulation of the basic 
Monte Carlo simulation routine is included below. 

7. 1 Basic Monte Carlo Simulation Routine 

The logic employed in the Monte Carlo simulation routine is shown in Figure 
7-1, To make use of this routine, the cloud-cover statistics (unconditional and con- 
ditional) are organized as cumulative probabilities in ascending order of cloud cover, 
such that the cumulative probability for cloud cover category 5 (total sky cover) is 
- 1. Mission itera'f'i number. Q, and pass number, n, are initialized. A random 
number, RAN, uniformly distributed in the interval O to 1 is generated and the 
first draw is made from the unconditional cumulative statistics by finding the cloua 
cover category having the cumulative probability interval containing the random 
number RAN. If the cloud cover category selected, G(n), is number 1 (i.e., clear 
skies), 100% observational coverage is assumed to have been achieved and the 
mission is completed. If the cloud cover category is not 1, then a percentage of 
the coverage, B, based on the cloud cover category, is assigned and the pass 
number is incremented by 1, A new cloud cover category is selected based on a 
second draw of random number RAN. This procedure progresses until the existing 
coverage is equal to 100% or that the number of passes equals a preassigned value 
N, This basic simulation scheme is quite flexible and has been successfully imple- 
mented by Greaves et al (1971), and others. Modifications to this basic scheme 
were made by Brown (1970) to provide a more realistic technique to increment the 
area of coverage after successive passes. 

7,2 Simulation to Include Cloud and Cloud Layer Statistics 

The suggested modifications to the basic simulation logic to include cloud 
layer and cloud type statistics are shown in Figure 7-2. The steps in this simula- 
tion scheme are identical to those of the basic scheme up to the point where the 
cloud cover category on the first draw is not number 1 and a percentage coverage, 

B, has to be assigned. At this point, a branch to ®*A** in the logic is made to 
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Figure 7-2 Modified Monte Carlo Simulation Routine 
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determine the %’crtical distiibutloua oi the cloud coverage by first determining the 
number of cloud layers and then the cloud type in each of the layers. 

In branch "A”, the statistics of cloud layers corresponding to the cloud 
cover category selected G(n) are first transformed into cumulative probabilities 
in ascending order of the number of layers. (Actually the order itself is not sig- 
nificant, ) A random number, RAN, is again drawn to determine the number of 
cloud layers. Again, RAN is uniformly distributed between 0 and 1 with the 
selection process for the cloud layers being identical to that for the cloud cover 
category. Having established the number of cloud layers, the process for the 
selection of the cloud type corresponding to each of the layers is initiated. As 
before, the statistics for the cloud types, UCT, corresponding to the given cloud 
cover category G(n) is first transformed to cumulative probabilities, CCT, in 
ascending order of the 1 to 9 cloud-type code used in this study. (Again, the 
order of the cloud types is not critical to the simulation scheme. ) The selection 
of the cloud type corresponding to the first layer proceeds through the same 
process of drawing a random number, RAN, uniformly distributed between 0 and 
1. The cloud type, CT, thus selected is then tabulated and its parameters such 
as height of base, thickness, droplet descriptors, can also be tabulated. 

At this point in the logic, the cloud-type selection is considered over if the 
number of cloud layers, NL, equals 1. If not a new cumulative probability for 
the cloud-type statistics is formed by subtracting from the total sample the fre- 
quency of occurrence of the cloud type, CT, just selected, thus in effect making 
the probability of selecting this cloud type in subsequent iteration zero. This 
step is suggested to ensure that two or more layers of the same cloud type are 
not generated in the simulation. The cloud layer number, NL, is then decreased 
by 1 and the process of cloud-type selection is repeated. The iteration continues 
until NL = 1 at which point the simulation returns to the original logic. 

7, 3 Discussion 

It is obvious from the above that the scheme is in its barest outline. A 
number of problems in the procedure have yet to be solved to make the simulation 
practical. These would include the question of spatial distributions of the cloud 
types and layers within the percentage of sky cover as specified by the cloud cover 
category and the logic for the specification of mission success or failure given the 
spectral dependence of cloud attenuation. 
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The basic data bank can and should be improved especially in the cloud-type 
distributions for those regions for which no taped data were conveniently available. 
The statistics of cloud parameters should also be updated on the basis of more 
data. This is especially necessary in the case of the ice-crystal and mixed clouds. 
With the data format adopted in this study, updating of the data bank to incorporate 
new data is a simple process. 
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APPENDIX A 


ANALYSIS OF THE EFFECTS OF RESOLUTION ON CLOUD 
COVER ESTIMATES 

The effects of sensor resolution size on cloud cover estimates derived 
from satellite photography were discussed in a study by Shenk and Salomonson 
(1971)^^\By means of a series of ingenious experiments, they derived a set of 
curves whereby estimated values of cloud cover can be obtained given the true 
cloud cover and a nondimensional parameter, R, defined as 

_ Areal Cloud Size 

Areal Resolution Element Size 

We have applied their results to an investigation of the effects of the parameter, R, 
on the conditional statistics of cloud cover derived by Sherr et al (1968)^^^ and 
extended by Greaves et al (1971)^^ ^ 

In this current effort, the emphasis was on estimating the true cloud cover 
from values of cloud cover derived from imperfect measurements; i. e. , AVCS- 
type photographs obtained by the Nimbus satellites. Satellite observations were 
the only source of data used by the referenced investigators to deduce the spatial 
conditional cloud cover statistics. In the past, these statistics were assumed to 
be perfect; i. e. , derived from observations made by a sensor with infinite spatial 
resolution, or at least with a resolution matching that of a surface observer. The 
resolution of the AVCS-type data is of the order of 1/2 mile. 


(1) Shenk, W. E. and V. V. Salomonson, 1971: A Simulation Study Exploring 
The Effects of Sensor Spatial Resolution on Estimates of Cloud Cover 
from Satellites , NASA T^ch. Note. Goddard Spaceflight Center. 
Greenbelt, Maryland. 

Sherr, P. E. , A. H. Glaser, J. C. Barnes and J. H. Willand, 1968: 

Worldwide Cloud Cover Distributions for Use in Computer Simulations, 
Final Report Contract NAS 8-2l040, Allied Research Associates , Inc . 

^ ^ Greaves, J. R. , D, B. Spiegler and J. H. Willand, 1971: Development 

of a Global Cloud Model for Simulating Earth- Viewing Space Missions, 
Final Report Contract NAS 8- 258 1 2^ Allfed Research Associates, 

Inc. 
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The results obtained by Shenk and Salomonson suggest that with the finite resolu- 
tion, there is a possibility that the cloud cover derived from the satellite data may 
in fact over-estimate the true cloud cover. For a given sensor resolution, the 
magnitude of this error is a fimction of the cloud size. As the cloud size decreases, 
the magnitude of this error increases. Figure A-2-1 shows this dependence. The 
curves arc derived from the basic data given by Shenk and Salomonson. The dif- 
ference lies in the fact that in this set of curves, the estimated cloud cover and 
the ratio, R, are considered to be the known parameters. 

It is clear that with this set of curves, the true cloud cover corresponding 
to each of the conditional cloud cover statistics can be derived if the following are 
known: 

1) Sensor Resolution 

2) Cloud Size 

While the former is known, the latter is not known in every instance. In fact, the 
best one can do is to obtain an estimate of the upper bound to the error by consider- 
ing the error obtained when the cloud size is represented by fair-weather cumulus 
which, as a cloud group, have the smallest horizontal dimensions. For this reason, 
a literature search was performed to obtain representative statistics of cumulus 
sizes. In this regard, a number of investigations were identified and in most of 
these the data were based on analysis of U-2 photographs. Presumably, the 
resolution of the camera systems are much better than the AVCS and are in fact 
assumed to be infinite. For this study, the data obtained by Blackmer and 
Serebreny (1962)^^^ were used to derive the characteristic dimensions of fair- 
weather cumulus. These data show that the widths of cumulus are characteristical- 
ly'^1 n. mi. Using this value and a sensor resolution of 2 n. mi. corresponding to 
the AVCS» the curve in Figure A -2 -2 was obtained which shows that at least for 
fair-weather cumulus cloud cover conditions, the conditional statistics derived 
from the AVCS data grossly under -estimate the true cloud conditions. It is im- 
portant to recognize that the data in Figure A -2 -2 represent an extreme condition. 
No set of spatial conditionals were derived using data containing fair-weather 
cumulus alone. For stratiform clouds it can be expected that R*^co and the error 
in the conditional statistics may be insignificant. 


^'^^Blackmer, R. H. and S. M. Serebreny, 1962; Dimensions and Distributions 
of Cumulus Clouds as Shown by U-2 Photographs, Report under Contract 
AF19(604)-7312, AFCRL-62-609, Stanford Research Institute. 
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Figure A-1-2 Example of Extreme Error in Cloud Cover Estimates 
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APPENDIX B 


SUBROUTINE AND CLOUD MODEL DEVELOPMENT 
The purpose of this Appendix is twofold: (1) To provide the listing 

of a simple subroutine for the computation of cloud-type and cloud-layer 
probabilities from the statistics compiled in the study and tabulated on 
cards and in Anpendix C; (2) to provide a set of more self-consistent 
models of liquid water clouds. 

1. SUBROUTINE CTPROB 

This is a general-purpose subroutine in FORTRAN IV, developed to 
compute useful cloud-type and cloud-layer probabilities using the tabulated 
cloud statistics. The subroutine is based on the simple arithmetic proba- 
bility computations discussed in Section 4 of this report. Specifically, 
the subroutine may be used to compute any or all of the following: 

a. The probability of occurrence of any of the nine cloud types 
or three cloud layer categories for any of the cloud-cover categories 
(except for the category of clear skies). 

b. The probability of occurrence of any of the nine cloud types 
and three cloud- layer categories independent of the cloud-cover category. 

c. The unconditional probability of occurrence of low (Cu, Sc, St), 

middle (As, Ac), and high (Ci, Cs) clouds, or combinations of any two or 
more cloud types. 

d. Unconditional probability of occurrence of each of the five 
cloud cover categories. 

The subroutine CTPROB requires four arguments: N, NF, P, CLRSKY, upon entry. 

Each of these is clearly identified and defined in the comment statements 
contained in the listing shown in the next two pages. 
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LISTING FOR SUBROUTINE CTPROB 

-Stf> IN 5 C Tt>rtr5-t W". NF". P . CL RSK Y t ' 

D PENSION NF(4,1? ) ,P(12> »lTEMPtI2) 




P ® D SE 

T.ALCtr.iTES P^:^34PILITIES OF CLDOD TYPES 4ND CL3UD L\Y = RS THF 

FRE3UI^-IES :3NT6!NE!) IN ARRAY NF, 


JSftGF 

- -CALu :TPR33tNtN=.P»CLRSKY> 

DISCRl’TION -OF PARAMETERS 

N -D^TION CD3E- SET N EQUAL TO 2.3«a,OR S T3 CALCJLATE )JD — 

TrP= -A^D- t-AYER PROftABlLlTI ES GIVEN ^LOUO“ COVER CAYE^ORY'V-. 

Spt'n equal TO 6 TO CALCULATE CLOUD TfPE AND LAYER 
P<33ABILITIE S INDEPENDENT OF CLFAR SKIES. 

S-T N =Q11AL TO 7 WHEN IT IS OESIR^D TO CALCULATE THE T3T AL 
j^j^^vj^-ITIlNAL probability OF OCCURRENCE QF CO'ASIJFD 

C.rn TYPES. THE COMBI NET CLOUD TYPES OF INTEREST 

— -"J-SYFT TN THF P^ ARHAY UPON ENTST, 

S=T N PQUAL TO 9 TO CALCULATE THE UNCONDITIONAL 

prtjability of OCCURRENCE of r . loud cover CATEGIRTFS, 

NF (-< t juc mE I: 3 OF OCC'.?PPENr=S OF 9 CLOUD TYPES 4N|T R CLOJO 
t AYERS FOR A CLOUD COVER CATEGORIES. 

P 'TNITAIMS THE RPSULTING PROBABILITIES GIVEN ^PFJ_ON N IN THE 

F3t-. OWI NG NANNEft . '■ ~ 

N=2.’.4.0I 5 - thru 9J=CL0UD TYPE PPOO ABIL IT I ES . 

odF THRO 12 I=CLOOO -layer DROBAR^ItlTIES. 

M = (S P SAME AS ABOVE. 

.^ = 7 Pm=TOTAL UNCONDITIONAL PROBABILITY F)E 

COMBINED TYPES. _ 

Y=9 - »-tt -T+tRU Sr-tJNCtM^OTTTONAtr-FROBA^tttT lES OF THE ' 

■5 CLOUD COVER CATEGORIES. 

CLRSSY the frequency OF CLEAR SKY. 


IFIN.GT.Bl S3 TO 11 
DO 10 1=1,12 


30 TO IS 

11 DO 12 1^1,12 

12 ITEMPCl )=0 
n DO 14 i=l.4 

DO 14 1=1,12 

- 1^1-TEMRN-H-s^tTENFf ! LfiNF FU. l »- 

15 TOTO=0.0 

on- 16 1=10,12 - 

16 T3TD=TOTO';FLOATnTFMPnn 
IF<N.GT.6)T0TD=TDTDCCLRSKV 

TFtN.'EQvB-lGO— TO— AG 

IF( N.E3.7I GO TO 20 

tot*o.o 

00 17 1=1,9 
RT I )=0.0 
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LISTING FOR SUBROUTINE CTPROB (contd) 


17 TOT = ij , , iLi'J (I» 

D ir 18 1 = 1,9 

IF( TOT.3T.O.O)® (I }*FL0AT( I TEMP( 1 ) »/TOT 
Ifl 

00 19 1 = 10,12 

I )^Q -.ir 

1 F( TDTD*3T,0.0»P« I )=F LCATdTEMP (I) ) /TOTD 

• 19 CJMTHUE — 

RETURN 

i-j(i) 

p a >=0.0 

21 - J=^,9 

IF( t» ( J) ,LT,1,0>GC TO 22 
< -P ( J ) 

21 ITEMP( I > = ITEMP( I) filTEMP(K) 



A sample run was made using the subroutine to ensure that it would 
execute properly and to illustrate its use with actual data. The data 
used were the cloud- type and cloud-layer statistics compiled for Month 1, 
time LST=1 for Region 19. These data are shown in Table 1, in which the 
frequencies of occurrence are presented in the format shown and discussed 
in Section 4 of the report. 

TABLE 1 

EXAMPLE OF CLOUD TYPE AND CLOUD LAYER FREQUENCIES 


Month 1 Region 19 


LST=1 

Cu 

Sc 

St 

Cb 

As 

Ac 

Ci 

Cs 

Ns 

1 

2 

3 

69 

2 

0 

3 

0 

0 

1 

5 

1 

4 

0 

10 

2 

0 


3 

0 

3 

0 

0 

0 

4 

1 

5 

0 

7 

3 

0 


4 

0 

1 

6 

0 

0 

1 

8 

1 

8 

1 

14 

5 

0 


5 

0 

28 

38 

0 

3 

11 

0 

13 

1 

51 

20 

1 



The results of the computation of probabilities, using CiPROB, for this 
set of input statistical values are shown in Table 2. 


2. FURTHER INVESTIGATIONS OF MODELS OF WATER CLOUDS 

It was noted in Section 6 of the Report that the drop size distribution 
of a liquid water cloud can be analytically expressed as: 


N(r) = a r^ 

where 

a 5 b are parameters specified for a given cloud type. 


( 1 ) 
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TABLE 2 

EXAMPLE OF CLOUD-TYPE AND CLOUD-LAYER PROBABILITIES COMPUTED USING CTPROB 
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According to Borovikov et al (1961) and others, these two parameters 

are related to the average drop radius, r , and to the liquid water content, 

W, as follows: 


a = 1.45 

(r)^ 


( 2 ) 


(r) 

It is clear from these expressions that self-consistent cloud models can 
be developed through the specifications of any two of the three parameters, 
N(r), W, and W appropriate to each of the cloud types. Since W and r are 
independent of the drop radius r, they are the more likely candidates as 
model descriptors. 

An examination of the recent literature on the microphysical properties 
of clouds shows that in developing empirical cloud- type models, the emphasis 
is almost always on the correct modeling of the drop size distribution with 
not much attention given to the modeling of the liquid water content. A 
very good example of this emphasis may be seen in the models developed 
by Cato et al (1965) . On the basis of analyses of a large number of drop 
size distribution measurements made by various experimentalists, they 
developed a series of cloud droplet size spectrum models, one for each 
of the major cloud types, for the purpose of theoretical investigations 
of scattering of laser light. The liquid water content of the cloud was 
considered not to be of importance. As a consequence, the models did not 
include specifications of this parameter. Despite this deficiency, the 
study of Cato et al provides a most complete set of cloud models based on 
actual measurements. We have, therefore, made use of these models as the 

starting point for the development of our analytical models. 
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Our approach towards the development of analytical cloud models was 
to preserve as much as possible the original empirical information contained 
in the models of Cato et al, and at the same time arrive at a set of para- 
meters such that Eq. (1) can be used. Values of r and W for each of the 
six types (no model of Ac was developed by Cato et al) of liquid water 
clouds were therefore computed using the ’’empirical” data of Cato et al 
and the following: 


r 

and 

W 

where 


r.N. 
1 11 


z 

i 


N. 


P i n (r.)- 

w 3 i ^ 


(3) 


(#/cm^ 


- ym) is the droplet number of density for droplet radius 

interval r. 

1 

is the density of liquid water 


with these values, the parameters a and b in Eq. (1) were computed by 
means of Eqs. (2) and (3). The resulting set of cloud model parameters 
is shown in Table 3. 

The computed values of a and b were used with Eq. (1) in sample 
calculations to determine how well the empirical distributions of N(r) 
are represented by the analytical expression. The results are shown in 
Table 4. 
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These data suggest that the analytic function fails at and near maximum 
values of N(r) when the distribution is sharply peaked as in the case of 
the stratus cloud model. Another feature of the models presented in Table 3 
which should be pointed out is that the liquid water content computed for 
the drop size distribution appears to be high. Figure 1, for example, 
shews the cumulative probability of occurrence of liquid water content 
measured in clouds. The curves are based on an analysis of the corresponding 
data given by Borovikov et al (1961). These data show that less than 20 % 
of the measurements show liquid water contents of greater than 0.3 gms/m^. 

Yet for three of the cloud models, Cb, As and Ns, the liquid water content 
is greater than 0.3 gms/m^. There is a possibility that the liquid water 
content model for As may be slightly too wet. 
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Figure 1 Cumulative Probability of Cloud Liquid Water Content 
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APPENDIX C 


TABULATION OF CLOUD- TYPE AND 

cloud-layeb statistics 
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